Rescaled density expansions and demixing in hard-sphere binary mixtures 
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■ The demixing transition of a binary fluid mixture of additive hard spheres is analyzed for different 
' size asymmetries by starting from the exact low-density expansion of the pressure. Already within 
' the second virial approximation the fluid separates into two phases of different composition with 

a lower consolute critical point. By successively incorporating the third, fourth, and fifth virial 
{2J[J' coefficients, the critical consolute point moves to higher values of the pressure and to lower values 

r-{ ' of the partial number fraction of the large spheres. When the exact low-density expansion of the 

pressure is rescaled to higher densities as in the Percus-Yevick theory, by adding more exact virial 
coefficients a different qualitative movement of the critical consolute point in the phase diagram is 
OO ' found. It is argued that the Percus-Yevick factor appearing in many empirical equations of state for 

T-H , the mixture has a deep influence on the location of the critical consolute point, so that the resulting 

phase diagram for a prescribed equation has to be taken with caution. 

^ I Demixing is a common phase transition in fluid mixtures usually originated on the asymmetry of the interactions 

a between the different components in the mixture. In the case of binary additive hard-sphere mixtures, the existence of 
I ' demixing has been studied theoretically since decades, and the issue has been at times controversial. The importance 
of this problem resides in the fact that if fluid-fluid separation occurs, it would represent a neat example of an entropy- 
driven phase transition, i.e. a phase separation based only on the size asymmetry of the spheres. The solution of the 
^ , Ornstein-Zernike (OZ) equation with the Percus-Yevick (PY) closure for additive hard-sphere mixtures by Lebowitai 
and the ensuing analysis of the thermodynamic properties of a binary mixture by Lebowitz and Rowlinson^ played a 
major role in the discussion and in many of the later theoretical developments. The great advantage of the PY theory 
C is that it yields analytic expressions for the equation of state (taking both the virial and the compressibility routes) 
as well as for the contact values of the radial distribution functions. Its main drawbacks are the fact that it predicts a 
(~| singularity in the pressure when the total volume fraction occupied by the spheres is equal to one (which is physically 
O unattainable) and the thermodynamic consistency problem, i.e. the fact that the virial and the compressibility routes 
, O, . lead to different expressions for the equation of state. In any case, the immediate conclusion that follows from this 
theory is that it predicts no separation into two fluid phases. The same conclusion is reached if one considers the most 
T-H , popular (and reasonably accurate when compared to simulation data) equation of state proposed for hard-sphere 
J> mixtures, namely the Boubhk-Mansoori-Carnahan-Starling-Leland (BMCSL)^ equation, which can be obtained by 
an empirical interpolation of the virial and the compressibility equations in the PY theory. However, the fluid-fluid 

■ segregation found by Biben and Hansen^ using the OZ equation with the Rogers- Young closure for highly asymmetric 
' spheres indicated otherwise and renewed interest in demixing. 

■ One can analyze theoretically the demixing transition in hard-sphere binary mixtures taking different routes. Of 
particular concern to us here is the approach based on the determination of an equation of state for the mixture 

, beyond the PY theory. Since no exact result is available, various possibilities arise, including the solution of the OZ 
• equation with different closures. However, this would most likely produce no analytical results. That is the reason 
'[ behind the many different empirical equations of state that have been proposed in the literature. Apart from aiming 
at providing a reasonably accurate account of the available simulation data, they attempt to reproduce virial behavior 
and/or to comply with consistency conditions of the contact values of the radial distribution functions (for a recent 
^ , review of these consistency conditions see Ref. |^ and references therein). In any event, all of them yield the exact 

second and third virial coefficients but at the same time inherit the singularity in the PY theory. 
Q , Coussaert and Baus^ have recently proposed an equation of state for a binary hard-sphere mixture with improved 
virial behavior that predicts a fluid-fluid transition at very high pressures (metastable with respect to freezing). On 
a different vein, Regnaut et al.^ have examined the connection between empirical expressions for the contact values 
of the radial distribution functions (which amounts to stating the corresponding equation of state for the mixture) 
5h ' and demixing. In this paper the demixing transition is analyzed by starting from the exact low-density expansion 
of the pressure. Our flndings are then compared with the fluid-fluid separation resulting from different empirical 
proposals for the equation of state. It is found that qualitative differences appear which confirm that demixing in 
binary additive hard-sphere mixtures is still an open question which deserves further investigation. 

We consider a binary fluid mixture oi N = Ni + N2 additive hard spheres of diameters txi and a2 (cri > 172). The 
thermodynamic properties of the mixture can be described in terms of the number density p = N/V, with V the 
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volume, the partial number fraction of the big spheres x = Ni/N, and the parameter 7 = (72/cri which measures the 
size asymmetry. For particles of equal mass, the Helmholtz free-energy per particle / = f{p,x,"f) reads: 

f3f = ln(/9A3) -l + a;lnx+(l-x)ln(l-x) 

00 ^ 

+ J2 -Bn+l{x,l)p\ (1) 
n—l 

where (3 = l/ksT, with T the absolute temperature and ks Boltzmann's constant, plays only the role of a scale factor, 
A is the thermal de Broglie wavelength of the particles, and Bn+i{x, 7) are the virial coefficients of the mixture. The 
pressure p = p{p, x, 7) is found from the thermodynamic relation p = p^(df /dp) yielding: 



(3p = p 



00 



l+> B„+i(x,7K 



71 = 1 



(2) 



By a Legendrc transformation we have moreover eliminated from (|5J| the number density in favor of the pressure and 
analyzed the phase separation using the Gibbs free-energy per particle g = g{p, 7) = / + p/p. 
For fixed 7, the lower critical consolute point, pc and Xc, is found from the convexity conditions: 
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Since g = xp,i -I- (1 — x)fi2, where /ii = fii {p, x, 7) and p2 = ^J■2{p^ x, 7) are the chemical potentials of the components, 
and from the Gibbs-Duhem relation x{dpi/dx) + (1 — x){dfj,2/dx) = 0, by also fixing the pressure {p > Pc), the partial 
number fractions of the big spheres at coexistence, x' and x" , are obtained from the coexistence conditions: 
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expressing the equality of the chemical potentials in both phases. 

We have truncated the density expansion in Eq. (j^J at fifth order since analytical expressions are knownS*^ for 
-82(2^,7) and 53(2:, 7), while Bn{x^^) and B^{x,'^) have been evaluated numericallyii2*ii*i24i^ Demixing has been 
investigated for small 7-values, say 7 = 0.05,0.1,0.2,0.3 and 0.4. For 7 = 0.1 we have also included a recent 
numerical evaluatior^ of Bq{x,j). 

We have first kept the density expansion up to second order, where the Gibbs free-energy can be obtained ana- 
lytically, and found that demixing occurs within this simple approximation. In Fig. 1 we plot the resulting phase 
diagram for different 7-values. We first note that the coexistence curves are restricted to a small region in the phase 
diagram since, for each 7, the pressure may only be increased until the total volume fraction occupied by the spheres, 
77 = 7r[x -I- (1 — x)j^]paf/6, reaches for the dense phase the highest possible value. Since these highest volume frac- 
tions are not known, all the curves have been drawn up to the limiting (unphysical) value rj — 1. It is seen that by 
decreasing the size asymmetry of the spheres the critical pressure first decreases and the coexistence region widens. 
After the critical pressure reaches a minimum (7 ~ 0.2) the trend is reversed. In contrast, the critical composition 
Xc always increases with 7. Note that when going for 7 = 0.3 to 7 = 0.4 the critical volume fraction rjc changes 
from 0.505 to 0.819, i.e. the coexistence region has shrunk substantially and this explains why by further increasing 
7 demixing disappears. We point out that this unexpected qualitative behavior is also encountered by keeping the 
virial expansion up to /cth-order (fc — 3, 4, 5) and similar to the phase diagram found in Ref. 0. 

An important feature for all the reported size asymmetries is that by successively incorporating the third, fourth, 
and fifth (sixth) virial coefficients, the critical pressure Pc increases whereas the critical partial number fraction Xc 
decreases, as shown in Table 1. 

The last finding seems to be in contradiction with the results by Coussaert and Baus who rescaled the density 
expansion in Eq. by assuming that the behavior of the pressure for high densities is the same as the one found in 
the FY theory, i.e. 
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FIG. 1: Phase diagram in the (a;,p*)-plane, with p* = f^po^i, of a binary additive hard-sphere mixture for different size 
asymmetries as obtained from the second virial approximation. The full dots denote the critical consolute points while the lines 
correspond to the coexistence curves. 
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Be 


7 


Pc Xc rjc 


Pc Xc rjc 


Pc Xc Tjc 


P*c Xc rjc 


Pc Xc Tjc 


0.05 
0.1 
0.2 
0.3 
0.4 


23.23 0.0304 0.276 
19.55 0.0368 0.274 
17.33 0.0561 0.339 
20.18 0.0869 0.505 
29.65 0.1300 0.819 


50.84 0.0294 0.428 
38.68 0.0365 0.394 
28.48 0.0525 0.381 
27.58 0.0769 0.447 
33.40 0.1134 0.580 


95.18 0.0254 0.530 
67.16 0.0333 0.476 
44.60 0.0498 0.432 
39.56 0.0727 0.463 
43.56 0.1066 0.548 


154.4 0.0223 0.598 
102.6 0.0308 0.534 
64.00 0.0489 0.476 
54.00 0.0724 0.492 
56.46 0.1054 0.554 


144.1 0.0289 0.576 



TABLE I: Critical constants p*, Xc, and rjc for different 7- values as obtained from the low-density expansion of the pressure by 
keeping n — 1 exact virial coefficients (indicated as Bn). 
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where the upper limit k in the sum indicates that Eq. Q guarantees, through C„+i(a;,7), that /3p has fc -I- 1 exact 
virial coefficients at low densities. Coussaert and Baus found that, when going from fc = 3 to fc = 4, the critical 
consolute point moves to lower values of pc and to higher values of Xc- Since the PY factor, 1/(1 — 77)^, strongly 
increases the pressure with respect to the pressure at low densities, the amplitude of the changes induced by adding 
one more exact virial coefRcient in Eq. ((SJ led the authors to conclude that the movement of the critical consolute 
point casts some doubt on the convergence of the procedure. Indeed, for a given k in Eq. JSj) the PY factor yields 
Bn{x,"f) ^ {n > k -\-2), which are otherwise zero in the exact density expansion up to (fc -I- l)t/i-order. Since 
the virial coefficients Bn{x,'y) (n > fc -f 2) in Eq. Q depend on x through the exact virial coefficients Bnix,^) 
{n < -I- 1), it can be concluded that the different additional dependences on x for fc = 3 and fc = 4 are responsible 
for shifting the critical consolute point to lower (higher) values of pc {xc)- Using a recent numerical evaluation of the 
sixth virial coefficient-^ for 7 = 0.1, we have confirmed the lack of convergence of this procedure and found that the 
critical consolute point moves from p* = 5377 and Xc = 0.0008 (fc = 3) to p* = 2565 and Xc = 0.0088 (fc = 4) and to 
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eCS-I 


eCS-II 
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P* Xc r]c 


Pc Xc rjc 


P* Xc rjc 


Pc Xc rjc 


0.05 
0.1 
0.2 
0.3 
0.4 


27.68 0.1120 0.589 
38.05 0.1114 0.624 
74.77 0.0244 0.389 
121.5 0.0263 0.471 
337.3 0.0263 0.620 


29.68 0.1001 0.595 

42.97 0.0968 0.636 

77.98 0.0197 0.366 

148.1 0.0211 0.478 

751.2 0.0179 0.681 


3599 0.0093 0.822 
1307 0.0203 0.757 
653.4 0.0537 0.725 
581.9 0.0998 0.738 
663.4 0.1532 0.766 


1096 0.0004 0.204 
832.0 0.0008 0.290 



TABLE II: The same as in Table I as obtained from the equations of state proposed by Hamad— (H-K and H-CS) and the two 
extended Carnahan-Starhng equations (eCS-I and eCS-II) of Santos et alt— 



p* = 1572 and x* = 0.0020 {k = 5). 

By the same argument, it is important to emphasize that the PY factor appearing in a great variety of empirical 
equations of state for a binary hard-sphere fluid mixture, arbitrarily introduces, as compared to the exact density 
expansion, additional x-dependences through the virial coefficients Bn{x,j) {n > 4) if, as a matter of fact, the 
empirical equations only preserve the first two exact virial coefficients. The additional x-dependences are also linked 
to the particular form of the empirical equation of state. This would explain why very similar equations of state do 
predict no demixing or demixing, and why, whenever it exists, the location of the critical consolute point strongly 
depends on the proposed empirical equation. In order to illustrate this point, we have analyzed four proposals for the 
equation of state that can be cast in the form of Eq. lO with either fc = 3 or A; = 4 and which predict demixingii^ 
They are the proposals by HamactlSi (which we will label as H-K and H-CS depending on whether the Kolafa or the 
Carnahan-Starling equation for the monocomponent ffuid is used as input) and the two extended Carnahan-Starling 
equations (eCS-I and eCS-II) of Santos et al?J^ All these equations have been introduced as alternatives to the BMCSL 
equation of state. It turns out that these four equations lead to the first two exact virial coefficients but differ in the 
predictions for Bn{x,j) (n > 4). Nevertheless they all give good agreement with the available simulation results for 
the compressibility factoriiSiiSiSSiSi The scatter in the values for the critical constants shown in Table II is so evident 
that there is no indication as to whether one should prefer one equation over the others, if any. For example, in the 
proposals by Hamad p* is always an increasing function of 7. In the cases of the eCS-I and eCS-II, p* does not follow 
this trend and can be up to two orders of magnitude higher. Moreover, the eCS-II does not predict demixing for 
7 > 0.2. 

In conclusion, we want to stress our main result, namely that the location of the critical consolute point of a 
binary hard-sphere fluid mixture, and hence the determination of the complete phase diagram (including the stability 
of demixing relative to freezing) have to be taken with caution when using empirical rescaled density expansions 
behaving in a PY-type manner. A deeper assessment is precluded at this stage and must await the determination of 
yet higher virial coefficients. 

After this paper was originally submitted, and due to a doubt that arose from the numerical value of the sixth 
virial coefficient as reported in Ref. IT3, we became aware through private correspondence with A. J. Masters of a 
complementary paper by Vlasov and Masters. There, the authors provide the most recent numerical values of the 
virial coefficients (up to the sixth) for various values of 7. For 7 = 0.1, they also compute the seventh virial coefficient 
and the critical consolute points considering different levels of truncation of the density expansion. Their critical 
constants are in good agreement with the data in Table I. 
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